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1. Introduction 

The theoretical proposal that Majorana modes exist in semiconductor devices [H EJ |3] 
and their subsequent detection in InSb wires [U El |6] has opened up a new subfield 
of research on nanostructure properties (see Refs. [3 El El [10] for reviews). It was 
originally proposed by Majorana that a massless elementary particle, called a Majorana 
particle, could exist with the peculiar property of being its own antiparticle. Similarly, 
a Majorana mode of a semiconductor nanowire is a zero-energy state that remains 
invariant after charge conjugation. These states are quasiparticle excitations localized 
on the tips of a finite but long enough wire and they are well separated from the rest of 
the spectrum of eigenstates by an energy gap. 

The existence of Majorana modes in a semiconductor wire requires the presence of 
the following physical ingredients: a) Zeeman coupling between spin and magnetic field, 
b) Rashba spin-orbit interaction and c) superconductivity [HI [121 02] ■ The latter can 
be induced by proximity with a superconductor material and it introduces the concept 
of electron- hole symmetry [T3]. The Rashba spin-orbit interaction is a relativistic 
effect originating in the quantum well asymmetry in the perpendicular direction to 
the nanostructure plane. In the present context this interaction introduces chirality 
by connecting the state of motion with spin. The Zeeman coupling in semiconductors 
like InAs and InSb is quite large even for relatively low magnetic fields due to the 
large g factors of these materials. It breaks Kramers degeneracy since the system is 
no longer time reversal invariant. In this paper we call Majorana nanowire (MNW) a 
semiconductor nanowire with all three physical effects a), b) and c) mentioned above. 

The transport properties in the presence of localized zero modes have been 
investigated for a normal-superconductor interface [151 UHl [T7]. It has been shown that 
both for resonant tunneling and for transmission through a quantum point contact a 
conductance quantization at half-integer multiples of Ae 2 /h is obtained when a zero 
mode is present at the interface. In this case the superconductor is called topological. 
In this manuscript we address a related although different geometry, the N/MNW/N 
structure where N refers to normal contacts in which the pairing and Rashba interactions 
vanish. We show that the existence of a zero mode is characterized by a unitary Andreev 
reflection, giving a linear conductance of the N/MNW/N structure equal to e 2 /h. We 
also find that by increasing the Zeeman coupling parallel to the wire a pronounced 
dip in the linear conductance appears due to the mixing between channels induced by 
the Rashba interaction. Similarly to the N/superconductor case |17] . the differential 
conductance has a peak at zero bias in presence of the zero mode, that evolves to a dip 
when the zero mode is absent. It is also worth stressing that the role of disorder of the 
N/MNW/N system has been studied in Ref. [T5] . 

Our main contribution in this work is the formalism of the coupled channel model 
(CCM) for the Bogliubov-deGennes (BdG) Hamiltonian. This formalism can be viewed 
as an alternative to the methods based on matching of plane waves, or on tight-binding 
chains [191 [20]. It is particularly adapted to the description of spatially smooth potentials 
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and it gives insights on the role of the different physical mechanisms by means of 
the channel-channel couplings. We present numerical solutions of the CCM equations 
for a representative case of a 2D MNW based on InAs. In support of our physical 
interpretations, we also present a simplified effective model allowing a fully analytical 
solution. 



2. The physical system 

The N/MNW/N system is modelled as a 2D channel of transverse dimension L y and 
with a central region of length L with superconducting and Rashba interactions. These 
interactions vary smoothly in the longitudinal direction taking constant values Ao and 
«o in the MNW, and zero in the asymptotic regions of the leads. In addition, potential 
barriers separate the central MNW from the leads. A sketch of the system and of the 
x-dependent functions is shown in Fig. [TJ The Hamiltonian reads 

UbcLG = (h - h)t z + A(x)r x + A B a ■ h 

H (PxO-y - PyO- x ) T z H — G y T z , (1) 

with 

2 2 

ho = + Pa + Va(x) + V c (y) . (2) 
2m 2m* 

The x-dependence of the pairing A(x), of the Rashba coupling a(x) and of the double 
barrier Vdb{x) is modeled by smooth Fermi-like functions with a small diffusivity d; for 
instance 

The transverse confinement potential V c (y) is taken simply as an infinite square well 
with zero potential at the bottom. The chemical potential explicitly appearing in the 
BdG theory is represented by parameter \i in Eq. ((T]); while a and r are, in a usual 
notation, the vectors of Pauli matrices acting in spin and isospin (or particle-hole) 
spaces, respectively. 

A distinctive characteristic of our model is the continuity of the system parameters 
with respect to the longitudinal coordinate x, as sketched in Fig. [Tb. This would allow 
us to investigate, for instance, the dependence on the diffusivity d of the transition. In 
this work, however, we will assume rather steep transitions of the system parameters. 
The coherent quasiparticle transport is described by the BdG equation 

U BdG ^ = E^> , (4) 

where E and \& are the quasiparticle energy and wave function, respectively. The latter 
depends on the position in space (x, y) as well as on the spin and isospin variables 
(VctjVt), where r\ =t, I represents a generic discrete variable with only two possible 
values, 

* => Vfay^rir) . (5) 
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Figure 1. Sketch of the physical system. Panel a) displays our choice of coordinates 
while panel b) shows the longitudinal variation of the Hamiltonian parameters. 

Notice, finally, that n is assumed to lie in the xy plane. An out-of-plane component 
would require the addition of orbital magnetic effects not considered in this work. 

3. The coupled channel model 

We present in this section the description of transport in terms of channel amplitudes or 
wave functions obeying a set of coupled differential equations. This description can be 
viewed as an alternative to the matching of bulk solutions often used in the literature. 
The CCM is well suited to the problem of a spatially continuous Hamiltonian posed 
in the preceding section. Before discussing the CCM equations, however, we need to 
consider the asymptotic solutions (x — > ±oo) as they are actually defining the channels 
themselves. 

3.1. Asymptotic solutions 

Since the pairing and Rashba intensities vanish asymptotically, the BdG Hamiltonian 
greatly simplifies in those regions, 

lim UbcLG = (h - h)t z + A B a ■ n . (6) 

In this limit the eigenstates are spinors pointing in the direction of n and z for spin and 
isospin. Introducing the quantum numbers s c = ±1 and s T = ±1 they read, respectively, 

Xsa ~ V2\ s a e^ J ' Xst ~ 2 ^ 1 - s T 

where ip is the azimuthal angle of h. The spatial dependence of the asymptotic 
eigenstates is also analytical, a plane wave in x and a square well eigenfunction in y, 
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4> n {y), with n = 1, 2, . . .. Summarizing, a channel is specified by the quantum numbers 
(ns a s T ) and its wave function reads 

y nSaSr = e ikx <p n {y)xsAv^XsAvr) • (8) 

The propagating or evanescent character of each channel is found when determining 
its wavenumber k = k nSaST . From Eq. (j3J) the asymptotic BdG energy is 



E =[^t^ + Zn-fi)s T + A B s a , (9) 



where 

ix 2 n 2 h 2 



2m*L 2 y 



(10) 



Inverting Eq. (J3|) it is 



/ 2 /77- * 

k nSaST = J — — (Es T - e n + /i - A B s a s T ) . (11) 

The channel wavenumber k nSiyST from Eq. (TTTT) is either real or purely imaginary. These 
two cases clearly correspond to propagating and evanescent channels, respectively. In 
conclusion, for electrons (s T = 1) and holes (s T = —1) the condition for propagating 
mode of spin s a in direction of h and with transverse state n is 

Es T - E n + [I - A B s a s T > . (12) 
3.2. The CCM equations 

Assume the following expansion of the full wave function, valid not only in the 
asymptotic leads but for any arbitrary position, 

^(x,y,Va,Vr) = ^ns a s T ( X ) ( l ) n(y)Xs tT (V<T)Xs T (Vr) , (13) 

where ip nSaST (x) is a ID function we call the channel amplitude. Obviously, the channel 
amplitudes in the asymptotic regions are VWst( x ) ^ ex P {^h nSaST x) : i.e., propagating 
or evanescent waves to the right or left directions, depending on the sign of the exponent. 
The equations fulfilled by the channel amplitudes can be obtained substituting the wave 
function, Eq. ( fT3l) . in the BdG equation, Eq. (BJ, and projecting on a specific channel, 

£ 1^ dy <f> n (y)x* s M)xt(r)r) x [H BdG * = EV] . (14) 

The sets of transverse wave functions {<p n }, {Xs a } and {xs T } fulfills proper 
orthonormality relations. After some straightforward algebra, Eq. ( |T4"1) leads to 



( V 2 \ 1 

( 7^ + V db (x) + e n - fi j s T + A B s a + s a s T smip— {p x ,a(x)} - E 

+s a s T cosip—{p x ,a(x)}i/j ns - aST (x) + A(x) i/j nSaS - T (x) 



-SaSr^- £ (n\p y \n') cos ip ip n i Sc ,s T (%) i sin ip ip n 's^s T {%) =0, (15) 

n'(^n) 
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where we have introduced the usual anticommutator notation, {p x ,a(x)} = p x a(x) + 
a(x)p x and the bar over an index denotes its opposite value. The set of Eqs. ffT5l) 
is already a first version of our desired CCM equations. There are three types of 
contributions to Eq. (fl5l) : a) the background terms of channel (ns a s T ) are given by 
the first line, b) the second line contains the coupling terms with channels of the same 
n but with opposite spin s c or isospin s T to that of the background channel, c) finally, 
the third line shows the coupling with channels of a different n, the same isospin and 
arbitrary spin. 

The physical role played by the different Hamiltonian contributions are clearly 
seen in Eq. (IT51) . As expected, the superconducting pairing A(x) couples electron and 
hole channels. The two Rashba terms have a markedly different effect regarding the 
n quantum number. The a(x)p x is diagonal in n, while the a(x)p y is mixing channels 
with different n's with the selection rules imposed by the square well matrix element 
(n\p y \n'). The relevance of the field orientation is also appreciated from Eq. (JT51) . For 
instance, if the field is along y {ip = tc/2) the mixing of (ns a s T ) and (ns~o-s T ) vanishes. 

We end this section mentioning a useful transformation of Eq. ffl5l) that eliminates 
the linear terms in p x of the background problem. Let us define the transformed channel 
amplitude 

i>ns a sM = e^ sin ^* (3;) Vw T (x) , (16) 

where we introduced the dimensionless function 

m* 



K. R {x) = I dx'a(x') . (17) 



o 



Substituting Eq. (fIS|) into Eq. (fl5|) we find 



+ 



+ V db (x) + e n - [X - ^ a ( x ) sin Vj s r + A B s a - E 
s T cos pe 2is ° sin v Kr{x) (i (x) 2 sin <p + s a ^ {p x , a (x) } 
A(x) i> ns<yS - T (x) - s a s T ^j^- ^ ( n \Py\ n ') cospijj n ' SrTST (x) 



- ismpe 2 ^^^ ij> n > 9 - aST (xj\ = , (18) 

The set of Eqs. (TlBl) is very similar to (fl5|) . with two important differences: a) 
the background channel terms have a new contribution quadratic in a(x) which is 
spin- independent, while the contribution linear in p x is effectively eliminated from this 
channel, b) the position-dependent phase of the transformation given in Eq. ffl6l) appears 
explicitly in the coupling with (nSo-s r ) and (n'so-s T ). 

3.3. The QTBM 

We have solved the set of Eqs. ( fT5i) using the quantum-transmitting-boundary 
formulation of the scattering problem. The reader is addressed to Refs. [22], [23] for 
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details on the QTBM. Here we just mention for the sake of completeness the basic 
underlying ideas. Using a ID grid Eq. ( jl5j) can be discretized with finite-difference 
formulas for the derivatives. In the asymptotic regions of the leads we impose the 
analytical solutions of the channel amplitudes 

lj) nf) s (x) = 0® e iSiS r kns a s T {x-Xi) _|_ e -iSiS T k„ a<7ST (x-Xi) ^ ^g-j 

where a^ SrT s r and b$ Sa s T are the usual incident and reflected amplitudes in lead i. In Eq. 
f|T9|) we have introduced the lead sign s iy equal to +1 and —1 for the left {i = 1) and 
right {i = 2) leads, respectively, as well as the position of each lead boundary Xi. We 
have also taken into account the reversed direction of propagation for electrons and holes 
with the s T sign. Notice that from Eq. ( IT9|) the outgoing coefficient bnl a s T is expressed 
in terms of the channel amplitude at the lead boundary, bnl aST = ipns a s r ( x i) ~ a nlas T - 
Substituting this explicit expression of bnl aST back in Eq. f fT9|) we obtain 

*Pns„s T (x) ~ e- ts ^ k "^ sAx - x ^tp ns<yST (xi) = 2ism( Si s T k nSaST (x - Xi))a® aBr . (20) 

The QTBM closed system of linear equations is defined as follows: a) for a grid 
point x such that x\ < x < X2 we impose the discretized version of Eq. f|T5|) . b) for a 
grid point having x < X\ or x > £2 we impose Eq. ( 1201) . The resulting linear system has 
as many equations as grid points and it depends only on the set of input coefficients 
{anl a s T }- It is highly sparse and can be numerically solved in an efficient way. 

The matrix of transmissions from mode ns a s T of lead i to mode n's' a s' T of lead i' is 
given by 



• 1 •/ / / / , 

t(ins a s T *>,., n ., TI - r (i) 



inSrrS T ) - 



ns a s T dns 



(21) 



where the subscript oim, standing for only incident mode, refers to the fact that all 
incident amplitudes vanish except the one explicitly appearing in the denominator of 
Eq. fT2~T|) . For use in the next section, we define a reduced matrix of transmission 
probabilities where we only discriminate lead and particle type, 

P°t= E Wns'J T ^in S(T s T )\ 2 . (22) 



4. Transport in the BdG framework 

The description of transport through MNW's can be done with the formalism of 
transport through normal/superconductor/normal structures. We follow, specifically, 
the formulation by Lambert et al. [21] for mesoscopic superconductors. For our two- 
terminal structure, labelled as % = 1, 2 for left and right contacts, the current in terminal 
i reads 

h = f°° E o ( - J t(E)) dE , (23) 

where E is the BdG quasiparticle energy. In Eq. ( 123|) . J"{E) and J^{E) are, respectively, 
the in- going and out-going fluxes in lead i of type a. 
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The essential ingredients we need to specify in order to use Eq. (1231) are the 
quasiparticle energy distributions f?(E), the number of propagating modes mf(E), 



and the matrix of quantum transmissions P°j (E). The latter two are obtained from 
the CCM, Eqs. (|T2l) and (|22|) . respectively. The quasiparticle distributions are assumed 
to be given by Fermi functions 

f?(E)= [i + e (£— vy/XT]-i ^ (24) 

where the i-th reservoir chemical potential has been defined as /ij = fi + eVi and kT is 
the thermal energy. With these inputs the fluxes in Eq. ( 1231) read 

J?(E) = ^m?(E)f?(E) , (25) 

JT(E) = ^PS P (E)ff(E). ( 26 ) 

jp 

This formalism fulfills two basic physical conditions: a) vanishing of current for 
zero bias and b) equality of current in both leads. Indeed, for zero bias all distributions 
are identical ft(E) = f(E) and then the sum rule on quantum transmissions, 

J2 P ?(E) = mUE), (27) 

ensures that in-going and out-going fluxes exactly cancel each other. The second 
condition, Ji + J 2 = 0, is more subtle; following Lambert [21] we interpret that it 
actually determines the MNW chemical potential /i, relative to /ii and /i 2 . Notice that 
the potential bias between the two leads is V = V\ — V 2 and that the MNW chemical 
potential lies somewhere in the range between the two reservoir chemical potentials, 

min(/ii,// 2 ) < \i < max(/ii,/i 2 ) . (28) 

The following practical approach to the BdG transport problem is then suggested: 
1) given Hi and /i 2 , assume fi = (/ii + /^ 2 )/2 and solve the BdG-CCM equations for the 
set {mf,P°j }; 2) compute Ji + I 2 ; 3) vary the value of /1 and recompute {m",P° /3 } 
until Ji + I 2 = is fulfilled. Solving this selfconsistency loop might be a difficult task, 
however it is not needed when the problem is symmetric with respect to x inversion 
around the center of the MNW. In this case /1 = (/ii + /i 2 )/2 is already the solution 
giving Ji + J 2 = since the bias V has to be shared symmetrically, V$ = SiV/2, where 
Si = 1 and s 2 = — 1. Here we shall focus on the symmetric problem, leaving for a future 
work the analysis on the non symmetric case. 



4-1- Differential and linear conductances 

The differential conductance, defined generically as dl/dV, is one of the most relevant 
transport properties usually measured in experiments. At zero temperature, the above 
formalism yields a very simple expression of this quantity because, in this limit, the 
derivatives of the quasiparticle distribution functions with respect to the bias become 
Dirac deltas. Of course, this is true only in the symmetric case, when V = 2s, V{. 
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For T = we obtain 

% = ^ (PiY(^V) + pr 2 -(hev) + Ptrtev) + P^eV)) , (29) 

and, as discussed above, it is d^/dV = —dI\/dV. The expression of the linear 
conductance G can be obtained simply setting the bias to zero in Eq. (129]) . Using, 
in addition, the particle hole symmetry 

Pff (E) = P?f{-E) , (30) 

we find 

G = j(P+ + (0) + P 1 \-(0)) . (31) 

Equations (I29|) and ( |3T]) are the basic relations of this work. Notice that 
they contain two qualitatively different contributions to the conductance, a normal 
transmission, T = (P^ + + P 12 _ )/2, whereby quasiparticle type is conserved; and an 
Andreev reflection, Ra = (P^T + Al + )/2) with quasiparticle change. Anticipating 
a result to be discussed below, notice that Eqs. ( 1291) and ( 13T1) predict a remarkable 
phenomenon, a nonvanishing conductance in absence of transmission (T = 0) due 
solely to Andreev reflection. This occurs when the Majorana nanowire has a zero mode. 
In this case Andreev reflection is maximal for zero bias, while increasing the bias there 
is a reduction of Ra, i.e., a zero-bias anomaly appears in dl\jdV due to the zero mode. 



5. Results and discussion 



5.1. Physical and scaled values of the parameters 

The relative strengths of spin-orbit, pairing and Zeeman terms for a given transverse 
dimension L y are determined by the following scaled dimensionless ratios (scaling is 
indicated with an s superscript) 

«? ' = Zg-L. , (32) 

A" = , (33) 

^ = ^2 ■ (34) 

Notice that for a given set of physical values of ato, A and A# different values of the 
transverse dimension L y will actually correspond to different relative strengths through 
Eqs. ( I32ll3^|) . Increasing L y the scenario clearly evolves from weak to strong couplings. 

More specifically, we consider below physical parameters that could represent an 
InAs-based nanowire j2l], m* = 0.033m e , cto = 30 meVnm and a pairing gap of A = 0.3 



meV. We assume a realistic value of the wire transverse dimension, L y = 150 nm, for 

(s) (s) 

which the relative strengths of Rashba and pairing are then a ~ 2 and Aq fa 3. 
Fixing these two scaled parameters to these values we will study the dependence on the 
third one A^ } below. The conversion of the Zeeman coupling to a physical magnetic 
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Figure 2. (Colour online) Linear conductance in absence of Rashba mixing as a 
function of the scaled Zeeman energy. The contributions from Andreev reflection (Ra) 
and normal transmission (To) are also shown. The inset shows a blow-up of the data 
in a small region. The system parameters are given in Subsection 5.1. 

field is B = Ag/g^B, with g and /i^ the g factor and Bohr magneton, respectively. 
With our assumptions this conversion reads B = (1.7 Ag /g) T, in terms of the scaled 
Zeeman coupling. That is, B = 1 T would correspond to A^ = 10 for a g factor w 17. 

The distance L between barriers (Fig. [1]) is taken as L = 20L y = 3 jum, with barrier 
thickness of 150 nm and height Vo = 0.5 meV, while the spatial diffusivity is d = 15 
nm [see, e.g., Eq. (j3J)]. We also choose the chemical potential, defining our reference 
energy of the MNW, as \i — 0. Overall, we stress that the complete parameter set is 
representative of a typical experiment with an InAs-based 2D semiconductor wire. 

5.2. Linear conductance results 

In Figs. |2] to H] we display the linear conductance calculated from Eq. ( )3T|) as a function 
of the scaled Zeeman value for the set of parameters mentioned in the preceding 
subsection. Figures |2] and |3] correspond to magnetic field in parallel direction to the 
wire (x) while Fig. 4 is for transverse orientation (y). In Figure 2 we neglected 
the contribution of the Rashba mixing, i.e., of the terms containing a(x)p y in Eq. 
( 1X5]) . Notice that in this situation the linear conductance displays an almost perfect 
quantization in e 2 /h steps. Increasing Ag, small deviations in the form of very narrow 
spikes can be seen at the beginning of the second and third plateaus. In the first two 
steps all the conductance is due to Andreev reflection since the normal transmission is 
negligible. Perfect Andreev reflection is a signal of the existence of a zero mode of the 
closed system, i.e., a Majorana fermion bound at the interface between the MNW and 
the normal contacts. This zero mode yields a perfectly quantized conductance in absence 
of transmission due solely to Andreev reflection, i.e., G = (e 2 /K)Ra. The decrease of 
Andreev reflection for A^ > 50 in Fig. [2] can be attributed to the finite size effect that 
removes the Majorana modes from perfect zero energy, in agreement with the analysis 
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Figure 3. (Colour online) Same as Fig. [2] including Rashba mixing. 




of Ref. [26J for the closed system. This decrease in Ra is accompanied by an increase 

in T , keeping the value of G close to integer multiples of e 2 /h, except at the transition 
between steps. 

Figure |3] displays the linear conductance for the same system of Fig. EJ but now 
including the full Rashba interaction. A conspicuous difference with Fig. [2] is that the 
conductance deviates from the simple staircase behaviour, with a broad conductance 
dip appearing at the end of the first plateau. This dip is due to a magnetic instability 
precluding the formation of two zero modes due to a repulsion between modes induced 
by Rashba mixing [26]. The effect of this mechanism on the linear conductance is 
remarkable, with the prediction of a reduced conductance due to a large reduction 
of Andreev reflection. This anomalous behaviour of the conductance at the end of the 
conductance plateau also appears in higher plateaus, as seen in Fig. [3] for the second and 
third plateaus. Notice, however, that the finite size effect mentioned above transforms 
the conductance dips of the higher plateaus in a strongly oscillating behaviour. We have 
checked that the formation of the conductance dips due to the instability of multiple zero 
modes is even more robust with higher values of the pairing gap and Rashba strengths, 
and that it is also robust against variations of the barriers between the normal contacts 
and the MNW. 

Figure H] shows the evolution of the linear conductance with for field along the 
transverse direction y. For this orientation of the field the physics changes completely, 
since now it is the Andreev reflection that vanishes and the conductance is due to the 
normal transmission. Only small peaks in Ra can be seen in the transition between 
plateaus. The vanishing of Ra is due to the absence of zero modes of the closed MNW 
for magnetic fields along y [26]. A similar orientation anisotropy has been seen in 
experiments with cylindric InSb nanowires [U E]. No conductance dips are observed 
in Fig. H] but there are many spikes due to resonant transmission through the double- 
barrier potential V^,{x). The separation between spikes is very small due to the dense 
distribution of quasibound states for such a long system L = 3 /im. From this point of 
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Figure 4. (Colour online) Same as Fig. [3] for magnetic field along y. 

view, it is still more remarkable the fact that for field along x the presence of a zero 
mode washes the spike oscillations and yields a consistent maximal conductance in some 
regimes. 

5.3. Nonlinear conductance 

The nonlinear conductance obtained with Eq. (1291) is shown in Fig. |5] as a function of 
the applied bias. We have taken some selected values of A^ from Fig. [31 corresponding 
to vanishing bias, and explored the variation with V. As in the preceding subsection 
we define a scaled bias taking the transverse confinement as reference, i.e., = 
(em*Ly/h 2 )V. For A^ = 9 there is a narrow peak in dli/dV at zero bias. This 
zero bias anomaly is reflecting the existence of a zero mode in the MNW. Increasing 
the Zeeman coupling the peak broadens, becoming a flat distribution. For A^ = 19, 
corresponding to the conductance dip of Fig. El the zero bias peak changes to a zero 
bias minimum. The existence of a zero bias anomaly in the presence of zero modes has 
been discussed before in systems with a superconductor contact, experimentally in Refs. 
[U E] and theoretically in Refs. [151 EH EE] ■ Our results prove that a similar behaviour 
is to be expected in N/MNW/N structures. 

5.4- Density distributions 

The density distributions, defined as \4 , nSa s T (x)\ 2 , are shown in Fig. [6] for two values 
of A^. They correspond to the perfect Andreev reflection (A^ = 10) and to the 
conductance dip (A^ = 19) of Fig. El As expected, the upper panel shows that the 
incident unitary density couples with an edge mode of the MNW. The density profile 
localized at the edge and decaying towards the interior has exactly the same shape 
found in calculations of zero modes of closed MNWs |26j. Perfect Andreev reflection in 
this situation consists in the total reflection in the conjugate channel and, therefore, no 
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Figure 5. (Colour online) dI\/dV as a function of scaled bias for the same 
parameters of Fig. [3] and for different values of A# . 
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Figure 6. Density distributions \ip nStrST (x)\ 2 for two values of of the first plateau 
of Fig. |31 We assumed boundary conditions corresponding to incidence from the left 
side. For comparison, the position of the potential barriers is indicated by the shaded 
regions. 



quantum interference is observed in the left contact. As mentioned before, this occurs 
due to the presence of the zero mode in the MNW and it allows unit conductance 
without any transmission at all between left and right contacts. 

The lower panel of Fig. O shows a qualitatively different behaviour. The beating 
pattern in the left contact is indicating that full reflection occurs now in the same 
channel of incidence, with a strong interference between incident and reflected waves 
The density at the edge of the MNW is more irregular and extends farther towards the 
interior than in the upper panel. The physical interpretation is clear: for this Zeeman 
intensity the edge mode of the MNW lies at a nonzero energy, this causing normal 
reflection, as opposed to the Andreev reflection of the upper panel. 
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6. A model in second quantization 

To provide additional insight on the physics of transport through the MNW, in this 
section we consider an effective Hamiltonian in second quantization projected onto the 
Majorana subspace. The Hamiltonian consists of three parts, that is, 

Heff = He + T~Lm + Ht , (35) 

with 

Tic = £ a kci k c ak , (36a) 

a=L/R,k 
1 

U M = ^mVlVr > (366) 
Kt = Yl {Kk,p c lkVfi + V ak) pr}pc ak ) . (36 c) 

a,P,k 

Here Tic describes the normal leads, with c ak (c ak ) the Dirac fermion creation 
(annihilation) operator. Notice that the spin degree of freedom is neglected. This can be 
understood considering that we need to apply a large magnetic field to observe the edge 
Majoranas, so that only one kind of spin is effectively involved |25j. Hm characterizes 
the coupled Majorana states, with T]l/r Majorana fermion operators fulfilling r/j = rj}, 
rjf = 1 and with anticommutator relation {f]i,r]j} = 25y. The parameter Em denotes 
the coupling between the two Majoranas on opposite ends of the MNW and can be 
some complicated function of wire length, superconducting coherence length, applied 
magnetic field, Rashba coupling and superconducting gap. Em might be found by exact 
diagonalization of Hamiltonian §Q (26] . We will assume that it is known for the purpose 
of the present model. The last contribution, %t, corresponds to the tunnel Hamiltonian 
between normal leads and the Majoranas on opposite ends |16j . Below, the tunnel 
amplitude V^p is taken as V for a = (3 and zero for a ^ (3. 
The current is computed from 



la = ~W,n*]) = ~»jEE^(^) \ , (37) 

P kGa 



where n a = J^kea c ak c ak and Q^k ,p denotes the lesser component of the mixed Green's 
function defined as 

Gak,p(t,t') = -i(T K c ak (t)r)p(t')) . (38) 

Employing the equation of motion technique, after tedious algebra, the current becomes 

la = ~ 9f || de Tr [/ a ( e ) (Q^e) - G a n {e)) T a {e) + Qf(e)T a {e)] } , (39) 

where G v {£) is the Green's function for the MNW and T a (e) denotes the hybridization 
matrix given by 

T a .p^(E) = it Y V akj pV* k J(e - E ak ) = 5 a> p5p n T/2 . (40) 
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To complete the calculation we need to determine the MNW Green's functions. They 
read [16] 

G r v /a (e) = , , (41) 

" e-it- 2^ a {e) 

G<(e) = g r Je)?:<(e)g"(e) . (42) 



Here, 



e M 
-£m 



(43) 



and the self-energy matrices are given by 

=o&(e) = =F<E ^(e) + r Tl *,(-e)] , (44) 

7 

S<^(e) = 2i £ l/ 7 (e)lWe) + / 7 (-e)r 7; ^(- £ )] . (45) 

7 

Substituting Eqs. (jlD]) - (J13|) into Eq. f )39|) and using current conservation we obtain 

h = -I R = e T l de — ^l r 2 2(£2 9 +4 r + 5 f) , : ^ (/*(*) - /*(*)) • ( 4 6) 



M (62 + 4^+6^(^-2(^-4^)) 

From Eq. ( 1461) we note that at T = the linear conductance finally reads 



e 2 4r 2 

G - J 4^ ' (47) 



For zero energy Majoranas £m — and then Eq. (jJBJ) yields G = e 2 //i. This result 
nicely agrees with our interpretation of Fig. |3] attributing maximal conductance to the 
zero mode. 



7. Conclusions 



We have presented the formalism of transport in a N/MNW/N structure based on 
the coupled channel model. This formalism yields a transparent interpretation of the 
coupling between channels induced by the relevant physical mechanisms of the problem. 
Namely, the confinement, Zeeman, Rashba and superconducting interactions. We have 
considered a 2D structure and in-plane magnetic fields, although the formalism can be 
extended to consider more spatial dimensions and different geometries. 

The coupled-channel-model equations have been solved using the quantum- 
transmitting-boundary algorithm for a set of parameters representative of an InAs 
nanowire. The existence of a zero mode in the MNW is characterized by a perfect 
Andreev reflection, whereby an incident channel is totally reflected in its antiparticle 
conjugate one. For a single zero mode the linear conductance takes the maximal 
value e 2 /h due solely to Andreev reflection, without any quantum transmission from 
left to right contacts. For increasing values of the Zeeman coupling along the wire, 
a conspicuous dip in the linear conductance is predicted due to repulsion between 



Transport through Majorana nanowires 



16 



Majoranas. This repulsion originates in the Rashba mixing between channels. On 
the contrary, for Zeeman coupling along y the Andreev reflection vanishes, with the 
possible exception of a small region close to the transition between plateaus. When 
the zero mode is absent, the linear conductance has narrow spikes as a function of the 
Zeeman coupling. 

The differential conductance signals the presence of the zero mode with a peak at 
zero bias. The zero bias peak evolves to a dip when the MNW zero mode is absent. 
Finally, we have also discussed an effective model in second quantization confirming 
the physical interpretation in terms of Majorana modes. The coupled channel model 
presented in this work can be used to investigate other scenarios like, e.g., non-symmetric 
barriers or sequential MNW's. Work along these lines is in progress. 
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